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Executive  Summary 


This  publication  presents  a  new 
approach  developed  for  analyzing 
geophysical  acoustic  reflection  data. 
The  forward  problem  of  estimating 
source  to  receiver  travel  time,  and  the 
inverse  probiem  of  estimating  a  sound 
speed  versus  depth  relationship  are 
addressed  using  a  moments  approach. 

The  moments  approach  provides  a  simple 
tool  for  estimating  source  to  receiver 
reflection  path  travel  time  for  a 
laterally  homogeneous  medium  with  an 
arbitvary  sound  speed  versus  depth 
relationship.  The  moments  approach 
also  provides  a  useful  tool  for 
inverting  reflection  data  to  obtain  an 
estimate  of  the  sound  speed  versus 
depth  relationship.  A  nearly  closed 
form  technique  for  estimating  a  linear 
sound  speed  versus  depth  relationship 
is  presented.  Derivation  of  the 
moments  approach  and  numerical  examples 
are  included. 
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A  Moments  Approach  for 
Analyzing  Geophysical  Reflection  Data 


I.  Introduction 

Geophysical  reflection  data  can  be 
analyzed  to  estimate  compress ional 
sound  speed  as  a  function  of  depth. 

Utility  of  such  data  may  be  subtle  or 
quite  straightforward;  exact  details 
are  beyond  the  scope  of  this  document. 
Analysis  presented  in  this  report  is 
confined  to  a  laterally  homogeneous 
med i urn  wi th  non-sloping  reflecting 
boundar i es . 

Acoustic  reflection  data  are  typically 
acquired  using  a  measurement  system 
geometry  similar  to  that  shown  in 
Figure  1.  The  source  provides  a  series 
of  acoustic  energy  pulses  which  prop¬ 
agate  through  the  medium  of  interest 
to  an  array  of  acoustic  receivers.  The 
distinctive  feature  of  acoustic  reflec¬ 
tion  data  is  that  the  received  energy 
arrives  via  a  reflection  phenomenon. 

The  reflection  phenomenon  is  caused  by 
an  abrupt  change  in  acoustic  impedance 
at  layer  boundaries.  Example  reflection 
paths  are  illustrated  by  paths  SR^  t> 

SR3  SRi  2  0^  Figure  1,  where  SRj  j 
denotes  the  reflection  path  from  source  to 
receiver  i  from  reflecting  boundary  j. 

Analyzing  acoustic  reflection  data  to 
determine  sound  speed  as  a  function  of 
depth  coordinate  in  Figure  1) 
requires  two  important  tools.  The 
first  tool  is  a  technique  for  estima¬ 
ting  reflection  path  travel  time  from 
source  to  receiver  for  candidate  sound 
speed  versus  depth  relationships.  The 
second  tool  is  the  inverse  of  the 
first,  in  that  the  second  tool  is  a 
technique  for  processing  measured 
travel  time  data  to  estimate  the  sound 
speed  versus  depth  relationship.  This 
report  describes  a  moments  approach  for 
handling  the  forward  problem  (travel 
time  from  known  sound  speed  profile) 
and  inverse  problem  (sound  speed 
profile  from  known  travel  time  data). 


The  forward  problem  is  handled  in  a 
general  sense  by  considering  an 
aTbitvavy  sound  speed  versus  depth 
relationship.  Treatment  of  the  inverse 
problem  is  not  so  general,  but  is  dis¬ 
tinguished  in  that  a  nearly  closed  form 
technique  of  estimating  a  linear  sound 
speed  versus  depth  relationship  is  pre¬ 
sented.  The  following  sections  of  this 
report  define  the  "moments"  of  a  sound 
speed  versus  depth  relationship  and 
illustrate  utility  of  the  approach  for 
processing  geophysical  acoustic  reflec¬ 
tion  data. 

II.  Moments  of  a  Sound  Speed  Profile 

The  compr es s i ona I  sound  speed  profile 
will  be  denoted  by  where  ^  is 

the  depth  coordinate  (y  =  0  at  the 
source).  As  d  i  scus seaear I  i  er ,  the 
analysis  presented  in  this  report 
considers  the  case  where  sound  speed  is 
a  function  of  depth  only  (laterally 
homogeneous  medium). 

Moments  of  a  sound  speed  versus  depth 
relationship,  Vzi^)  >  *  particular 

measurement  configuration  are  defined 
with  the  aid  of  Figure  2.  The  i*^ 
moment*,  Mj ,  is  defined  as  follows. 


'Moments  of  a  sound  speed  vs.  depth 
relationship  should  not  be  confused 
with  moments  of  a  probability  density 
function  (statistics)  or  moments  of  a 
mass  distribution  (mechanics).  The 
particular  label  "moments"  was  chosen 
because  of  slight  similarities  to 
statistical  and  mechanical  definitions. 


ARRAY  OF  ACOUSTIC  RECEIVERS 


Figure  1 .  Measurement  System  Geometry 
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Figure  2.  Defining  Moments 
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+  V^’(D^  -  D,)  for  i  ;»  -1  . 


where  ^  =  depth  below  source, 

V2(*)  =  sound  speed  as  a 
function  of 

Dj^  =  *  coordinate  of  reflecting 
Doundary  and 

Df  =  ^  coordinate  of  acoustic 
receiver. 


where  £n(.)  denotes  natural  log. 

Additional  examples  are  given  in 
Appendix  A.  The  important  feature  to 
be  noted  is  that  moments  can  be  cal¬ 
culated  quite  easily  for  virtually  any 
Vz(jj)  (sound  soeed  vs.  depth). 

III.  Estimating  Source  to  Receiver  Reflection 
Path  Travel  Time  (Forward  Problem) 


For  the  case  where  source  and  receiver 
are  at  the  same  depth  (D,.  =  0)  ,  the 
moments  definition  reduces  to 


u 

M.  A2  / 


'^2 


Without  regard  to  their  utility,  it  is 
beneficial  to  notice  the  simplicity  of 
computing  moments  for  a  sound  speed 
profile  Vz(^).  As  an  example,  con¬ 
sider  the  source  and  receiver  geometry 
of  Figure  3  v/ith  a  Vz(j^)  consisting 
of  a  constant  gradient  (  \  intiv  sound 
soeed  versus  depth  relationship)  layer 
and  a  constant  sound  speed  layer.  From 
defining  equation  (1b),  flj  for  the 
case  displayed  in  Figure  3  becomes 


M,  =  2 


y  u, 

J  (Vo.g;)id;*2  /  V^d;.  (2) 

0  Di 


A.  Problem  Definition 

The  "forward  problem"  associated  with 
analyzing  reflection  data  is  defined 
as  estimating  the  source  to  receiver 
reflection  oath  travel  time  for  a  given 
measurement  geometry  and  Vz(*).  A 
sample  measurement  geometry  ms  shown  in 
Figure  3.  In  this  case  (Fig.  3),  the 
forward  problem  is  to  estimate  the 
source  to  receiver  reflection  path 
travel  time  as  a  function  of  source  to 
receiver  horizontal  "offset"  distance 
X.  The  travel  time  as  a  function  of  x, 
T(x),  is  cormonly  labeled  "moveout". 

Computing  T(x)  for  a  geometry  consist¬ 
ing  of  a  single  constant  sound  speed 
layer,  as  shown  in  Figure  4a,  is  quite 
simnie  and  results  in  the  hyperbolic 
"moveout"  T(x)  as  shov/n  below  and  in 
Figure  4b. 


T(x)  =  (r(0)  + 


Evaluating  equation  (2)  results  in  the 
f 0  I  lowing. 


2(0.  -  D,) 

* 


for  i  =  -1  and 


where  x  =  source  to  receiver 
hor i zon  t  a  I  offset, 

A 

Vg  =  sound  speed  in  layer 

(constant  in  this  case) 
and 

T(0)  =  Normal  incidence  reflec¬ 
tion  path  travel  time. 
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SOURCE  receiver 


Figure  4a.  Single  Constant  Sound  Speed  Layer 


X 

Figure  4b.  T(x)  for  Single  Constant  Sound  Speed  Layer 


5 


Com<-!U  t  i  n"  T(x)  for  a  general 
much  more  difficult.  The  most 
straightforward  approach  for  computing 
T(x)  for  a  general 

trace  technique,  which  approximates 
^z(j)  *  series  of  simple  functions 

ancT  then  incoroorates  Snell's  law  to 
determine  an  approximate  travel  time 
and  distance  as  a  function  of  sound  ray 
departure  angle.  Error  in  the 
approximation  can  be  made  small  at  the 
expense  of  computational  burden  by 
increasing  the  number  of  simple  func¬ 
tions  approximating  It  is 

precisely  the  computational  burden 
associated  with  ray  tracing  that 
motivated  develoonent  of  the  moments 
approach  for  estimating  T(x). 

B.  Moments  Approach  for  Estimating  T(x) 

The  moments  approach  for  estimating 
"moveout",  T(x),  anproximates  T(x)  by 
the  square  root  of  a  finite  order 
Maclaurin  series  about  the  normal 
incidence  reflection  path,  viz. 

T(x)  =  (C  +  Cx  +  Cx  ^  Lx  +  C.x'  — )'^  (3) 


discussed  later,  is  that  it  lends 
itself  more  readily  to  the  inverse 
problem. 

Clearly  the  challenge  associated  with 
implementing  equation  (4)  is  deter¬ 
mining  the  coefficients  Cj .  A  prin¬ 
cipal  feature  of  the  moments  approach 
is  that  it  allows  the  coefficients, 

Cj  ,  to  be  computed  quite  easily. 

Coefficient  Cq  can  be  related  to 
physical  conditions  by  evaluating 
equation  (4)  at  zero  offset  (x  =  0), 
i  *  c  • } 

c  =  T(x  -  0).  (5) 

Referring  to  the  general  configuration 
of  Figure  5  leads  to  the  follov/ing 
equation  for  T(x=0). 

■  /'  vf/T 

0  ^  D  ^  ' 


where  x  =  source  to  receiver 

horizontal  offset  and 


Cj  =  function  of  measurement 
geome  try  and  . 


Ry  syrrmetry,  it  is  easily  shown  (May 
and  Straley  Ml)  that  odd-ordered  coef¬ 
ficients  C-j  ,  Cg ,  etc.,  are  identically 
zero.  Deleting  zero  valued  coefficients 
and  squaring  equation  (3)  results  in 


T  (x)  =  C  +  C  X  +  C  X  +  C,  X  - .  (4) 


The  series  approximation  is  no  dif¬ 
ferent  from  that  of  Taner  and  Koehler 
12]  and  flay  and  Straley  |1  |.  The 
moments  approach  differs,  ho\/ever,  in 
that  it  addresses  the  ‘  ■>•1:1  ‘  V'l  rj  vz(t) 
as  opposed  to  layers  i-f  constant  soxind 
speed.  A  further  advantage  of  the 
moments  approach,  which  v/ill  be 


Comparing  equations  (6)  and  (la) 
reveals  that  T(x=0)  and  therefore  Cg 
can  be  expressed  identically  in  terms 
of  moments  by 


C  =  (fl_^)' 


(7) 


Differentiating  equation  (4)  with 
respect  to  x^  and  taking  the  I imi t  as 
X  goes  to  zero  results  in  the  following 
express  ion  for  C2- 


C  =  Lin, 

X  *0  dx' 


(8) 


Evaluating  dT^/dx^  in  terms  of 
physical  parameters  can  be  simplified 
by  applying  a  mathematical  identity  to 
obtain 

■dj.  .  T  _dj/dx 
dx  ^ 


(9) 


Figure  5.  General  Configuration 


Also,  dT/dx  can  be  expressed  by 


dT  _ 


dx  dx/dt)., 


(10) 


Referring  to  the  notation  defined  in 
Figure  5,  expressions  for  x  and  T  can 
be  derived  to  yield 


X  =  /  tan  ■i{^) 


(lla) 


sinu(y)  =  sinUj,  and 


do 


VAf)  cosOq 

v'^oT~cos0ijy 


(13a) 


{13b) 


Substituting  e'^uations  (13)  into  equa¬ 
tions  (12)  results  in  the  following 
expressions  for  dx/d0o  and  dT/d9o 


n 

dx_^  , 

I/'  f  Jiz'r)  J 

j"  tan  .(j)  df 

and 

V^fO) 

y  cos^o(^)  cos‘y(^) 

T  .  /■ 


f 


S'Ht) 


(Ub) 


where  6q  =  initial  departure  angle 
of  sound  ray  and 


and 


dX  "’""o  cos  0,^1 


do. 


V^'IO) 


■¥ 


cos 


dj 


I 


Vz(/) 


0^  cos  \A}) 


(14a) 


(14b) 


9(«)  =  sound  ray  angle  at  depth^ 
for  a  ray  departing  v/ith^ 
angle  Oq' 

Oi c farent i a t I nj  equations  (11a)  and 
(1  b)  \/ith  respect  to  Qq  results  in 


rinci 


D 

f  " 

Q 


. —  a* 

cos  (jt)  ' 


D  d<;(r) 

(it  +  J  — . -  d^  ( l?a) 


0. 


cos  I) 


(;) 


Substituting  equations  (14)  into 
equation  (10)  yields  the  following 
familiar  expression  (Rutherford  131) 


dT 

dx 


sin'.i 

v^Tci) 


(25) 


Subs t i tu t i n"  equations  (13)  into 
equation  (11a)  results  in  the  following 
equation  for  x  as  a  function  of  9q. 


i 


4 

•? 


% 

1 

i 


dT 
d  ' 


+ 


sin  'i^)d 

- ( - d  * 

V^(/)co3  -■{j! 

. . . d* 

Vz(/)cos  ' 


From  Snell's  law,  the  following 
relationships  can  be  obtained. 


(  12b) 


> 


sin 


/ 


v/^i 
(os’  r^') 


cos  (^1 


J 


If, 


An  cxnression  for  dlZ^rjx^  can  now 
be  obtained  by  substituting,  enuations 
(16),  (15),  and  (lib)  into  equation 
(9) ,  viz. 


« 


-  J  *diS 


A.1  interesting  by--iro'tuct  of  this 
result  is  that  the  Dix  14)  approach  for 
riet  erni  n  i  rig  flflS  sound  soeed  of  a  layer 
is  exact  (as  opposed  to  aijprox  inate) 
for  an  arbitrary  ,  providing 

that  "array  velocity^',  is  defined 

as  foil  ov/s  . 


To  evaluate  coefficient  C2  (equation 
(R)),  v/e  need  the  limit  of  dT^/dx^ 
as  x  aoproaches  zero.  Realizing  that 
Gq  approaching  zero  is  equivalent  to 
X  approaching  zero,  allows  equation 
(17)  to  be  evaluated  in  the  limit  to 
yield 


Lim 

X  >0  dx' 


Comparing  equation  (18)  v/i  th  defining 
equation  (la)  reveals  that  equation 
(18)  and  therefore  coefficient  C2 
(via  equation  (8))  can  be  expressed 
identically  in  terms  of  moments  by 


u  Ci  Lim 
''a  '  x-0 


Higher  order  coefficients  of  the 
t2(x)  polynomial  are  derived  in 
Appendix  For  brevity,  only  the 

results  are  presented  here. 


C  -  (M_,)' 


(20a) 


(20b) 


M  M 

_ I 

Mr' 


(20c) 


C.' 


(19) 


2 


M/  M 

-mT' 


-  M, 


M,  -  M,  M., 


(20d) 


It  is  easily  shown  (see  section  IV. 
0.5)  that  is  identically 

the  Rf'iS  sound  speed  for  the  normal 
incidence  (x=0)  reflection  oath  from 
source  to  receiver,  where  the  RfiS 
speed,  Vpji,|5,  is  defined  by 

Tn  j  V.(t)dt 
0  ' 

where  Tn  =  Normal  incidence  reflec¬ 
tion  path  t  r  ave I  t ime 
and 

Vj ( t )  =  Sound  speed  as 

f unc  t i on  of  t ime . 


'^RMS 


C.  - 

64M, 


9  M,  M,  *  24  M,  M5  M 


- 1 


-  4  M,  M.-'  -  5  M,  M.,  M; 


24  M.  M_ 

m; 


(20e) 


The  technique  presented  in  Appendix  0 
can  be  used  to  derive  higher  order 
coefficients  (C^q,  etc.).  Practical 
experience,  however,  has  indicated  that 
an  eighth  order  "moveout"  t2(x)  poly¬ 
nomial  is  more  than  adequate  to  obtain 
accuracy  comparable  to  a  practical  ray 
trace  a  I gor i t  hm. 
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The  important  feature  of  this  deriva¬ 
tion  is  that  the  finite  order  t2(x) 
polynomial  can  be  generated  easily 
using  the  moments  approach  for  vir¬ 
tually  any  of  interest. 

C.  Accuracy  of  the  Moments  Approach  for 
Estimating  T{x) 

Accuracy  of  the  moments  approach  for 
estimating  reflection  path  "moveout” 
was  investigated  by  comparing  T(x) 
computed  from  the  moments  approach, 
T|n|(x),  with  that  computed  from  a 
standard  ray  trace  algorithm,  Tf{(x). 
Experiments  were  conducted  for  a 
variety  of  sound  speed  profiles,  a  few 
of  which  are  presented  here. 

The  first  experiment  considers  a  single 
200  m  thick  layer  with  a  sound  speed 
gradient  of  1.5  sec"^.  Figure  6a 
displays  the  measurement  configuration 
and  sound  speed  profile.  Figure  6b 
displays  the  error  (difference  between 
ray  trace  and  moments  approach)  as  a 
function  of  source  to  receiver  hor- 
i,!ontal  offset  x.  Errors  for  the 
second  order  through  eighth  order 
moments  approach  are  presented.  Notice 
that  an  eighth  order  approximation 
estimates  T(x)  to  within  2  micro¬ 
seconds  . 

The  second  experiment  involves  a  more 
detailed  sound  speed  profile  as  shown 
in  Figure  7a.  Performance  of  the 
moments  approach  is  illustrated  in 
Figure  7b . 

The  third  experiment  is  designed  to 
show  ihvaat  L  /  ■  ty  and  corraut at  tonal 
s  nnnyi'  offered  by  the  moments  approach 
for  handling  a  more  general  A 

quadratic  sound  speed  versus  depth 
relationship,  as  shown  in  Figure  8a,  is 
considered.  In  order  to  implement  the 
ray  trace  algorithm,  it  was  necessary 
to  lyi  noxr'hi' ■'  the  quadratic 
a  set  of  N  linear  sound  speed  versus 
depth  relationships.  Accuracy  of  the 
ray  trace  algorithm  improves,  of 
course,  with  increasing  N.  Figure  8b 
displays  the  difference  (T(x)  esti¬ 
mates)  between  ray  trace  and  the  eighth 
order  moments  aproach  as  a  fijnetion  of 


N  for  a  source  to  receiver  offset  of 
9S0  m.  Notice  that  as  N  becomes  large, 
the  ray  trace  estimate,  T[^,  converges 
near  the  mome nts  approach  estimate. 
Computational  savings  offered  by  the 
moments  approach  over  ray  trace  for 
N=100  is  approximately  two  orders  of 
magn i t  ude . 

iV.  Estimating  Vz(^  from  Reflection  Data 
(Inverse  Problem) 

A.  Problem  Definition 

The  "inverse  problem"  of  estimating  a 
sound  speed  versus  depth  relationship 
from  acoustic  reflection  travel  time 
data  is  considerably  more  difficult 
than  the  forward  problem.  The  level  of 
difficulty  increases  with  generality  of 
the  solution.  As  an  example,  the  Dix 
|4j  approach  for  estimating  RMS  sound 
speed  of  a  layer  is  quite  simple, 
whereas  estimating  a  linear  sound  speed 
versus  depth  relationship  (Gibson, 
Odegard  and  Sutton  15])  can  be  orders 
of  magnitude  more  difficult.  The  Dix 
approach,  therefore,  features  simplic¬ 
ity  at  the  expense  of  only  estimating 
RfIS  sound  speed  (1  parameter)  of  the 
layer  of  interest.  The  technique  used 
by  Gibson  et  al.|5l  features  more 
generality  (linear  sound  speed  versus 
depth  relationship)  at  the  expense  of 
additional  computational  burden.  The 
particular  approach  used  by  Gibson  et 
al.  j5|  uses  a  ray  trace  technique 
coupled  with  a  nonlinear  estimation 
algorithm  to  determine  layer  character¬ 
istics  (thickness,  gradient,  etc.) 
necessary  to  matah  the  observed 
moveout,  T(x).  The  moments  approach 
inversion  technique  presented  in  this 
report  features  generality  (linear 
sound  speed  versus  depth  relationship) 
with  only  a  slight  increase  in 
computational  burden. 


*Actual  results,  using  available 
software,  indicated  a  400  to  1  time 
reduction  in  computational  burden 
using  the  moments  approach. 
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Figure  6a.  Measurement  Configuration  and  Sound  Speed 
Versus  Depth  Relationship 


Figure  6b.  Errors  in  T(x)  Computed  from  Moments  Approach 


DEPTH  BELOW  SOURCE  AND  RECEIVER  (METERS) 


|Tg(X)-T^  (X)|(SECONDS) 


B.  Moments  Approach  for  Estimating  a 
Linear  Sound  ^eed  Versus  Depth 
Relationship  (Single  Layer  Case) 


squared,  Performing 

this  step  results  in  a  polynomial 
expression  as  shown  below, 


The  moments  approach  offers  a  very 
simple  technique  for  estimating  a 
linear  sound  speed  versus  depth  rela¬ 
tionship  from  acoustic  reflection  data. 
Specifically,  the  approach  estimates 
the  following  parameters  describing  the 
I ayer  of  interest. 

o  Thickness  (H) 

o  Sound  speed  at  "top"  of  the  layer 

(Vt) 

o  Sound  speed  at  "bottom"  of  the  layer 
(Vr) 

o  Sound  speed  gradient,  dV^/dj^,  (g) 

0  RMS  sound  speed  (Vf^j) 


^  Co  ^  .  C„xO  ... 

where  T,„eas(*)  =  Measured  reflection 
path  travel  t  ime  as 
a  function  of  source 
to  rece  i  ver  hor i zon- 
tal  offset,  X. 

Many  techniques  (least  squares  etc.) 
exist  for  determining  the  coefficients, 
Cj .  Details  of  these  techniques  will 
not  be  described  here.  For  a  linear 
sound  speed  versus  depth  relationship, 
only  coefficients  Cg,  C2  and  C4 
are  needed. 


For  simplicity,  the  a  I  go r i t hm  wi I  I  be 
first  derived  for  a  single  layer  case. 
Also,  for  simplicity,  the  algorithm  is 
derived  for  the  case  where  source  and 
receiver  (or  an  array  of  receivers)  are 
at  the  s  ame  depth.  Me  a  s  u  r  erne  n  t 
geometry  and  V2(^)  for  the  single 
layer  case  are  shown  in  Figure  9.  The 
source  and  receiver  (or  array  of 
receivers)  is  used  to  measure  the 
reflection  path  travel  time  T(x)  where 
X  is  the  source  to  receiver  horizontal 
offset. 


(Compute  Moments' 


The  second  step  relates  moments  of  the 
sound  speed  profile  to  coefficients, 

Cj  ,  determined  in  step  1.  For  a 
linear  sound  speed  versus  depth  rela¬ 
tionship,  only  moments  M_i  ,  M^j  and 
M3  need  be  computed.  Equations  (20) 
can  be  solved  to  yield  the  following 
expressions  for  relating  moments  to 
polynomial  coefficients. 


H 


(Cr) 


1 

2 


The  procedure  for  estimating  layer 
parameters  from  measured  T(x)  (denoted 
by  T^eas(*))  very  much  the  reverse 
of  estimating  T(x)  from  known  layer 
parameters.  More  spec i f i c a  I  I y ,  the 
procedure  follows  in  steps  as  outlined 
be  I ow. 

Step  1;  Relate  moveout  polynomial 
coefficients,  Cj.  to 
measured  travel  time  data. 


Step  2:  Relate  moments,  Mj ,  to 
coef f i c i ents  Cj . 

Step  3:  Relate  layer  parameters,  H, 
Vj ,  etc.,  to  mome n t s . 


3.  S  t  ep 
Sol ut ion] 


(Closed  Form  Approximate 


The  third  step,  relating  layer  par am- 
e  t  e  r  s  to  mome  n  t  s ,  i s  s  omewh  at  more 
difficult.  Relating  moments  to  layer 
par  ame  t  e  r  s .  for  a  linear  sound  speed 
versus  depth  relationship,  using  equa¬ 
tion  (1b)  results  in  the  following. 


The  first  step  is  simply  "fitting"  a 
polynomial  to  the  measured  moveout 
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Figure  9,  Single  Layer  with  Linear  Sound  Speed  Versus  Depth  Relationshi 


M,  =  (Vg  -  V^) 


«3  =  (v;  -  v;) 


where  Vj  =  Sound  speed  at  "top"  of 
the  layer, 

Vg  =  Sound  speed  at  "bottom"  of 
^  the  1  ayer , 

g  =  Sound  speed  gradient  and 

iln(,)  denotes  natural  logarithm. 

Substituting  for  g  (g  =(V0  -  Vj)/H) 
where  H  is  the  layer  thickness,  results 
in  a  more  convenient  form. 


Unfortunately,  equation  (24)  involves 
the  unknown  parameters  Vj  and  Vg. 
Equation  (24)  will  be  used,  however,  in 
the  itepative  solution  presented  in  the 
next  section.  For  practical  cases,  n 
will  be  very  nearly  unity. 

Equations  (23),  (22b)  and  (22c)  can  be 
solved  in  closed  form  to  yield  the 
following  expressions  relating  layer 
parameter  estimates  (H,  etc.)  to 
mome  n  t  s . 


Vt  (n) 


Vg  (n)  = 


M,  +  Q 


M,  =  (V^  +  V3)H 


A  /Ms  M-j 
where  Q  =  Kl - 


(Vt  +  VJ  (VJ  +  Vp')H 


K  =  +  1  (See  Section  B.6, 
page  20 ) 


Inverting  equations  (22)  to  find  layer 
parameters  as  a  function  of  moments, 
involves  the  solution  of  nonlinear 
algebraic  equations.  Unfortunately,  no 
closed  form  solution  has  been  found.  A 
alosed  form  approximate  solution  can  be 
obtained,  however,  by  expanding  the 
natural  log  function  of  equation  (22a) 
about  Vj/Vg  =  1  (Ref.  [6|).  Per¬ 
forming  this  expansion,  and  introducing 
a  correction  factor,  r| ,  results  in  the 
following  representation  for  equation 
(22a). 


The  closed  form  approximate  solution 
assumes  n  is  identically  unity.  With 
this  assumption,  the  closed  form 
approximate  solution  becomes 


Ml  :  -Q 

2H 


Ml  t  Q 


The  correction  factor,  r|f  is  given  by 


where  Q  =  K  (M3  M_'|  -  and 

K  =  +  1  (See  Section  B.6, 
page  20) 

4.  Step  3  (Iterative  Solution) 


^"(vr)<''T"  V 


For  many  applications,  the  closed  form 
approximate  solution  of  equations  (25) 
provides  acceptable  results.  If 
greater  accuracy  is  desired,  the  ttera- 
tivc  solution  described  in  Figure  10 
can  be  used. 


The  "loop"  from  Block  to  Block  (?) 
of  Figure  10  is  repeated  until  esti¬ 
mates  for  Vg,  and  H  converge. 

As  an  example,  consider  a  layer  des¬ 
cribed  by  the  following  parameters. 

H  =  100  m 

Vj  =  1 500  m/sec 
Vg  =  1600  m/sec 

Table  1  illustrates  properties  of  the 
iterative  solution  by  tabulating  H, 

Vj,  and  Vg  as  a  function  of  the 
number  of  iterations.  Results  of  the 
closed  form  approximate  solution  are 
shown  as  zero  iterations. 

5.  Es  t  imat  i  ng  Gr  ad  i  ent  and  R^.IS  Sound 
Speed 

Once  Vj,  Vg  and  H  are  determined, 
an  estimate  of  sound  speed  gradient 
follows  directly  from  the  definition  of 
gradient  for  a  constant  gradient  layer, 
i  .e. , 


g  =  (Vg  -  v^)/H 


(26) 


The  RMS  sound  speed,  of  the 

layer  is  defined  as  the  root-mean- 
square  speed  of  a  sound  ray  travel¬ 
ling  at  normal  incidence,  i.e.. 


(t)  dt 


(27) 


where  Tp  =  Normal  incidence  (one 

way)  travel  time  from^=0 
tojf=D£  of  Figure  9  and 

Vj(t)  =  Speed  of  a  sound  ray 
t  r  ave I i ng  at  norma  I 
incidence  as  a  function 
of  time,  t. 

Rearranging  the  integration  of  equation 
(27)  allows  VR)»i)5  to  be  expressed  by 

''RMS  =  (f  f  d(t(y)))  (28) 

^  0 

The  function  Vy(t  (/))  is,  of  cour  se , 
identical  I  y  V2_{^) . 

Vy  (t(;))  =  V^C;)  (29) 

Also  the  function  t(*)  can  be  expressed 
by 


"  j 

0 

Differentiating  equation  (30)  with 
respect  to  *  and  solving  for  the  dif¬ 
ferential  d(t(^))  results  in 

d(t(^))  =  '  (31) 


Table  1.  Example  of  Closed  Form  Approximate  ond  Iterative  Solution  (Single  layer  Case) 


yv  Av 


Number 

of  1  terat ions 

H  (m) 

Vt  (m/sec) 

Vr  (m/sec) 

0 

(Closed  form 

100.017 

1491.994 

1607.468 

1 

approx,  so  1  ut  ion) 

99.994 

1502.958 

1597.222 

2 

100.002 

1499.048 

1600.892 

3 

99.999 

1 500.322 

1599.698 

4 

100.000 

1499.893 

1600.100 
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Substituting  equations  (31)  and  (29) 
into  equation  (28)  results  in  the 
following  expression  for  VRjy|5. 


Relating  equation  (32a)  to  moments 
yields 


which  is  identically  (€2)'^  • 

Therefore,  an  estimate  of  RMS  sound 
speed  of  the  layer  can  be  obtained 
inmediately  after  Step  1,  viz. 


where  Vj  =  Sound  speed  at  "top"  of 
layer  and 

Vg  =  Sound  speed  at  "bottom"  of 
I  ayer 

Notice  from  equations  (33)  that  if  Vi¬ 
and  Vg  are  interchanged  (g  becomes  -g) 
the  moments,  Mj ,  remain  unchanged. 
Consequently,  the  moveout,  T(x),  for 
reflection  path  data,  is  not  a  function 
of  the  polarity  of  the  sound  speed  gra¬ 
dient.  Therefore,  successful  "inver¬ 
sion"  depends  on  a  priori  knowledge  of 
the  polarity  of  sound  speed  gradients 
in  constant  gradient  layers.  In  view  of 
this,  the  rules  for  picking  K  of  equations 
(25)  and  Figure  10  bee ome 

K  =  •^1  For  sound  speed  i  nc  r  eas  i  ng  with 
increasing  depth  and 

K  =  -1  For  sound  speed  dec r ea s i ng  with 
increasing  depth. 


^DM<;  ■  !  C  )  ' 


where  C2  is  a  coefficient  in  the 
moveout  polynomial.  Equation  (32b)  is 
true  for  any  sound  speed  versus  depth 
relationship,  whereas  Equation  (26)  is 
true  only  for  a  constant  gradient  case. 

6.  Amb i gu i ty 

An  unfortunate  ambiguity  exists  when 
"inverting"  reflection  path  travel  time 
data  from  a  constant  gradient  layer. 

To  explain  the  ambiguity,  equations  for 
moments  of  a  constant  gradient  layer 
are  r  epe  a  ted  here. 


7.  Sunmary  (Single  Layer  Case) 

For  clarity,  the  procedure  for  esti¬ 
mating  a  linear  sound  speed  versus 
depth  relationship  from  reflection  path 
travel  time  data  is  sumnarized  here. 


Step  1 :  Fit  polynomial  to  measured 
moveout  data  Tff,gj5(x)  to  obtain  coe 
ficients,  Cj  ,  of  the  moveout  poly¬ 
nomial 


.  (x)  =  Cfl 


Cux''  +  Ce,x^ 


Step  2 :  Compute  moments  from  moveout 
polynomial  coefficients. 
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Step  3  (Closed  Form  Approximate 
Solution):  Compute  estimates  of  layer 

par  ame  t  e  r  s  f  r  om  mome  n  t  s . 


where  Q  =  K{l1i  M_ ,  -  M]')  ^ 

K  =  +t  for  sound  speed  i n c r e a s - 
i n g  with  increasing  depth 

K  =  -1  for  sound  speed  decreas- 
i ng  with  increasing  depth 

Step  3  (Iterative  Solution):  To 
improve  accuracy  of  Step  3,  the  itera¬ 
tive  solution  described  in  Figure  10 
may  be  used. 

Estimating  gradien t  and  Rf,1S  sound  speed 
is  performed  very  simply  by: 

9  =  (Vg  -  V^)/H 

Vrms  =  (C  ) 

C.  Moments  Approach  for  Estimating  a 
Linear  Sound  Speed  Versus  Depth 
Relationship  (Multilayer  Case) 

Extending  the  mome nts  approach  to 
estimate  a  linear  sound  speed  versus 
depth  relationship  for  the  multilayer 
case  is  reasonably  simple.  Basically, 
the  extension  involves  processing  tu'O 
measured  moveouts,  T{x),  to  determ. ne 
moments  for  the  layer  of  interest. 

Figure  11  illustrates  the  general  meas¬ 
urement  configuration  and  a  multilayer 
sound  speed  profile.  Moveout  data  for 
the  upper  reflecting  boundary,  Tu(xJ, 
and  lower  reflecting  boundary,  ( x )  , 
are  used  to  compute  moments,  Mj ,  of 
the  layer  of  interest.  The  moments, 

Mj  ,  are  then  used  to  estimate  layer 
parameters  by  incorporating  the  singte 
Ici'ijCV  procedure  derived  in  previous 
sect  ions. 


The  procedure  for  estimating  layer 
parameters  of  interest  follows  in  steps 
and  is  described  below. 

1 .  Step  1  (Determine  "Upper"  Moments) 

Observe  moveout,  Tu  (x),  of  the 
reflection  path  associated  with  the 
upper  reflecting  boundary  (see  Fig. 

11).  Compute  polynomial  coefficients, 
Cuj,  describing  measured  moveout  to 
satisfy 


Tu  (yl  =  Cu  I  +  Cu  X  +  Cu  X  ■  +  Cu  x’ 
meas 


Compute  "upper"  moments.  Mu j ,  from 
coefficients  Cu j .  As  in  the  single 
layer  case,  only  3  moments  are  required 
for  the  linear  sound  speed  versus  depth 
relationship. 


Mu  .  =  (Cu  ,)^ 


(34a) 


Mu- 


(34b) 


(34c) 


2 .  Step  2  (Determine  "Lower"  Moments) 
Observe  moveout,  I^.,,(x),  of  the 

meets ' 

reflection  path  associated  with  the 
lower  reflecting  boundary  (see  Fig. 
11).  Compute  polynomial  coefficients, 
Ck.j,  describing  measured  moveout  to 
satisfy 


Compute  "lower"  moments,  Micj  ,  from 
coef f i c  i  ent  s  CX j . 

Mi.,  =  (CD,,)''  (35a) 
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Figure  11.  Multilayer  Sound  Speed  Profile  (General  Model) 


MV,  =  (l  -  ^ 

\Ci2/  V  Ci. 


3.  Step  3  (Compute  Moments) 

Compute  moments,  Mj ,  for  the  layer  of 
interest  using  the  upper  and  lower 
mome  n  t  s . 


M.  =  Mjc^  -  Mu.  i  =  -1,1,3 


4.  Step  4  (Closed  Form  Approximate 
Solution] 


Relate  layer  parameters  of  interest  to 
moments  computed  from  step  3. 


(V-)'’ 


.  M,  -  Q 


V  =  ^  Q 

®  2H 


g  =  (Vg  -  v^}  /  H 


where 


Q  =  K  (M,  M.  I  -  M,-'!  ' 


K  =  +1  for  sound  speed  increas¬ 
ing  with  increasing  depth 


K  =  -1  for  sound  speed  decreas¬ 
ing  with  increasing  depth 


5.  Step  4  (Iterative  Solution] 

In  many  cases,  the  closed  form  approx¬ 
imate  solution  of  equations  (37)  is 
accurate  enough.  If  better  accuracy  is 
desired,  the  iterative  solution  of 
Figure  10  can  be  incorporated.  In 
either  case  (closed  form  approximate  or 
iterative)  estimates  of  sound  speed 
gradient,  g,  and  RMS  sound  speed, 

Vf^S>  given  by  equations  (38)  and 

(39). 

6.  Derivation  of  Step  3 

Steps  1,  2  and  4  are  essentially  iden¬ 
tical  to  their  single  layer  case  coun¬ 
terparts  and  therefore  will  not  be  dis¬ 
cussed  fur  ther . 

Step  3,  specifically  equation  (36), 
will  now  be  derived.  The  moments  of 
interest  describe  the  layer  of  interest 
and  are  therefore  defined  by 


M.  =  2 


where  Dy  and  Dj^  are  depths  of  the 
upper  and  lower  reflecting  boundaries 
(see  Fig.  11).  The  "upper"  moments, 
Mu j  ,  computed  in  step  1  describe 
for  .  The  "I  owe  r " 

moments,  Mj^j  ,  computed  in  step  2 
describe  V^i*)  for  O^^Dj^.  There¬ 
fore,  Muj  ano  M)lj  re^esent  the 
following  information. 


.  yu  . 


U 

r  f 


Subtracting  equation  (41a)  from  (41b) 
therefore  results  in  the  moment:  of 
interest,  i.e.. 
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.  -  Mu^ 


(42) 

7.  Multilayer  Ex  amp  I  e 

Figure  12  displays  an  example  measure¬ 
ment  system  geometry  and  multilayer 
sound  speed  profile.  The  layer  of 
interest  is  described  by 

H  =  60  m, 

Vj  =  1 560  m/ sec ,  and 
Vg  =  1 650  m/ s ec . 

Steps  1,  2  and  3  were  executed  to 
obtain  three  moments  (M_i  ,  M-|  , 

M3)  for  the  layer  of  interest.  Step 
4  was  executed  to  obtain  the  closed 
form  approximate  and  iterative  solution 
for  layer  parameters.  Results  of  step 
4  are  displayed  in  Table  2.  Results  of 
the  closed  form  approximate  solution 
are  displayed  as  0  iterations. 

8.  Features  of  the  Solution 

Notice  from  the  derivations  that  com¬ 
putations  for  Mj  depend  only  on  move- 
out  data  from  reflecting  boundaries 
directly  above  and  below  the  layer  of 
interest.  Therefore,  the  moments 
approach  for  estimating  a  linear  sound 
speed  versus  depth  relationship  is 
similar  to  the  Dix  (4|  approach  for 
estimating  Rf15  sound  speed.  The  sig¬ 
nificance  of  this  similarity  is  that 
only  errors  (measurement  or  otherwise) 
associated  with  two  sets  of  measurement 
data  influence  the  parameter  estimates 
for  the  layer  of  interest. 


Also  notice  that  while  a  linear  sound 
speed  versus  depth  relationship  V^l^) 
is  assumed  for  the  layer  of  interest, 
absolutely  no  oonstraining  assumption 
is  made  concerning  outside  the 

layer  of  interest. 

V.  Summary  and  Recommendations 

A  "moments"  approach  has  been  developed 
for  processing  geophysical  acoustic 
reflection  data.  The  major  feature  of 
the  moments  approach  is  a  tremendous 
reduction  in  computational  burden  over 
standard  ray  trace  techniques.  The 
"forward"  problem  of  computing  reflec¬ 
tion  path  travel  time  has  been  solved 
using  the  moments  approach  for  the  case 
of  an  arbitrary  sound  speed  versus 
depth  relationship.  The  "inverse" 
problem  of  estimating  a  sound  speed 
versus  depth  realtionship  from  measured 
reflection  data  has  also  been  addressed 
using  the  moments  approach.  A  closed 
form  approximate  solution  and  an 
iterative  solution  has  been  developed 
for  estimating  a  linear  sound  speed 
versus  depth  relationship  (constant 
gradient).  The  derivation  makes  no 
constraining  assumption  concerning  the 
sound  speed  versus  depth  relationship 
outside  the  layer  of  interest.  All 
derivations,  forward  and  inverse, 
assume  a  laterally  homogeneous  medium 
with  nonsloping  boundaries. 

Additional  efforts  should  be  directed 
toward  extending  the  moments  approach 
to  sloping  boundaries  and  perhaps 
laterally  nonhomogeneoiis  meiHiuns. 

Also,  an  extension  to  higher  order 
(higher  than  linear  sound  speed  versus 


Table  2  Example  of  Closed  Form  Approximate  and  Iterative  Solotian  (Multiloyer  Case) 


Numb  er  of  Iterations 

H  (m) 

Vt  (m/ sec ) 

Vr  (m/sec) 

0  (Closed  form 

60.008 

1 552.827 

1656.753 

approx,  solut  ion) 

I 

59.997 

1562.649 

1647.491 

2 

60.001 

1559.148 

1650.805 

3 

60.000 

1 560.288 

1649.728 

4 

60.000 

1559.904 

1650.090 
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depth  relationships)  inversion  tech¬ 
niques  should  be  investigated  along 
with  error  analyses  for  all  reflection 
data  inversion  techniques. 
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Appendix  A. 

Moments  for  Exemple  Sound  Spood  Profilos 


I.  Series  of  Constant  Sound  Speed  Layers 

A  sound  speed  versus  depth  relation¬ 
ship,  V2(x)i  consisting  of  N  constant 
sound  spefd  layers  is  shown  in  Figure 
A-1 .  From  defining  equation  (la)  of 
the  text,  moments  Mj  of  the  sound 
speed  profile  (for  the  given  measure¬ 
ment  configuration)  are  given  by 


For  the  case  where  receiver  depth  is 
same  as  source  depth  (D]=0)  equation 
(A. 2a)  reduces  to 


(A. 2b) 


Substituting  the  moments  of  equation 
(A. 2b)  into  equations  (20)  of  the  text 
results  in  moveout,  T(x),  polynomial 
coefficients  which  agree  (to  within 
obvious  typographical  errors)  with  that 
reported  by  Taner  and  Koehler  (Alj. 

il.  Series  of  Constant  Gradient  Layers 


(A.l) 


where  is  receiver  depth  and  D{v| 
is  depth  of  the  reflecting  boundary  of 
i nteres t .  Substituting  the  sound  speed 
profile,  V2(*),  illustrated  by  Figure 
A-1  into  equation  (A.l),  results  in  the 
f 0 1  lowing. 


A  sound  speed  profile,  (V2(*) ,  consist¬ 
ing  of  N  constant  gradient  ^\iyeri  is 
shown  in  Figure  A-2.  From  defining 
equation  (la)  of  the  text,  moments,  Mj  , 
of  the  sound  speed  profile  (for  the 
given  measurement  configuration)  are 
given  by 


^  Vj  if)  d;  *  /"n  vj  (,)  <1, 

0  0, 


(A. 3) 


N  D. 


°J-. 


Sound  speed,  Vj ,  being  constant  over 
the  specified  integration  intervals 
allows  Mj  to  be  expressed  by 


"i  -i:  ■  »j-,)  *  z  'joj  -  »j.,i 

j=l  j=2 


where  Di  is  receiver  depth  and 
is  depth  of  the  reflecting  boundary  of 
interest.  Adopting  the  notation  Vt j 
and  Vb j  are  velocities  at  the  "top" 
and  "bottom",  respectively,  of  layer  j, 
allows  the  moments  to  be  expressed  by 


(A. 4a) 


and 


Rearranging,  and  noticing  (Dj  -  Oj_ij 

is  identically  Hj ,  results  in  the  more 
convenient  form 


^  =  TTPTyE 

j=i 


_ J  - 


M,  =  2 


N 

E 

j=i 


^  - 


v;  H. 


{A.2a) 


(Vb’/‘  -  Vt’/') 
91  (1  +  1) 


for  i  t  -\  . 


(A.4b) 
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SOURCE 

* 

^  RECEIVER 


1 

Vz 

r^vt, 

LAYER  1 

Vvb, 

GRADIENT  =  g 

LAYER  2 

,v 

GRADIENT  =  g 

K. 

LAYER  3 

GRADIENT  =  g 

GRADIENT  =  9,^  LAYER  N 

\ 

REFLECTING  BOUNDARY 

^OF  INTEREST 

Figure  A-2.  Series  of  Constant  Gradient  Layers 
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For  the  case  where  receiver  depth  is  the 
same  as  source  depth  (D^=0)  equatiorts 
(A. 4)  reduce  to 

N  [Vb-I 


,  N  Vb’^  '  Vt’*‘ 

^  ^  rr^£ 

J=-  ^ 

III.  More  General  V2(j;) 

Con  luing  the  list  of  examples  is 
som  lat  pointless,  since  it  is 
do  .tful  that  an  example  wilt  fully 
m<;et  the  requir  erne  nts  of  an  actual 
application.  Montents  for  more  general 
sound  speed  versus  depth  relationship 
can  be  computed  directly  from  defining 
equations  (1)  of  the  text. 

REFERENCES 
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Appendix  B. 

Derivation  of  T^(x)  Poiynomiai 


,  D  Vj(^)  if  Vj{^)  d/ 

E.  .  =  f  - —  +  /  - —  (6.3) 

’’J  Q  cosJ"(^)  •'  cos^oC;?) 


I.  Relating  Coefficients  to  Derivatives 

The  "moveout"  polynomial  squared  is 
approximated!  by  a  finite  order  Taylor 
series  about  the  normal  incidence 
reflection  travel  path;  viz. 


T(x)  =  C,  +  Cx"  +  C.x"  +  C^x^  +  ... 


where  t  =  Depth  below  source  [see 
Fig.  B-11 

Vz(«)  =  Sound  speed  at  depth  « 

'  [  see  F i g.  B-1  ) 


Initial  sound  ray 
departure  angl e  | see 
Fig.  B-1 1 


Coefficients  of  the  expansion  can  be 
related  to  T^(x)  and  its  derivatives 
to  yield  the  following  equations. 


=  Lim  T'  (x) 


C  Lim 


f  Lin, 


1  d  r  dT' 


dx  I  dx 


J  cIT 

^  dx  dx  dx 


(B.2a) 


(B.2b) 


(B.2c) 


(B.2d) 


Q{m)  =  Sound  ray  angle  with 
'  respect  to  vertical  at 
depth  ^  for  an  initial 
departure  angl e  ©q 
I  see  Fig.  B-1 ) 

A 

Djj  =  ^  coord  i  na  t  e  of 

reflecting  boundary  of 
interest  |see  Fig.  B-l] 

A 

Dr  =  ^coordinate  of 

receiver  j  see  Fig.  B-1  | 

Differentiating  equation  (B.3)  with 
respect  to  ©q  and  applying  Snell's 
law  |Vz(0)  sin  O(^)  =  Vz(jr)  sin  ©q  1 
results  in  the  following  useful  rela¬ 
tionship. 


'i.j  ^  ^  .J* 


1  0 

^  N  ^ n;  :r~ji 


(B.2p) 


where  Q 


sin  cos 


II.  Notation 

Deriving  the  derivative  expression  of 
equations  (B.2)  in  terms  of  physical 
considerations  (measurement  geonretry 
and  sound  speed  versus  depth,  Vz(«)) 
IS  greatly  simplified  by  defining^the 
following  notation: 


Reflection  path  travel  time  and  x  direc 
tion  travel  distance  as  a  function  of 
sound  ray  Initial  departure  angle  *-'0 
can  be  expressed  by  | see  Fig.  B-1 | 
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SOURCE 


gure  B-1,  General  Configuration 


x(U(,) 

=  /■“’’  sinoM 

J  COS0(*) 

0  ' 

sina(* 
J  cosa(t 
Dr 

j-d^  (6.5) 

=  f  * _ ^ 

d« 

i  ^ 

r 

c5s0(^) 

(B.6) 

Applying  Snell’s  law  1V2(0)  sin  0(*) 
Vz(/)  sin  0Q  I  to  equation  (B . S  ) 
yields 


(B.7) 


receiver  x=0 )  travel  time  squared. 
Therefore,  Cq  can  be  expressed  as 
foil ows . 


Equation  (B.13)  can  also  be  expressed  in 
terms  of  moments  by 

C,  = 

where  moments,  Mj ,  are  defined  by 

D, 

'  v’(/)d^  (B.14) 


Using  the  notation  of  defining  equation 
(B.3),  equations  (B.6)  and  (B.7)  and 
their  derivatives  with  respect  to  Oq 
can  be  expressed  by 

\  1  t  (B.8) 


IV.  Determining  C2 

Coefficient  C2  ‘s  expressed  in  terms  of 
dT^/dx^in  equation  (B.2b).  The  deri¬ 
vative,  dT^/dx^,  can  be  expanded  (by 
identity)  to  yield 


sin  cos 

.)i  whore  0 

V^(0) 


dr 

dx 


T_dJ/^x 

X 


(B.15) 


sin 

V^'(O)  '• 


Substituting  equations  (B.8),  (B.10) 
and  (B.12)  into  equation  (B.15)  yields 
the  following  expression  for  dT^/dx^. 


cos 

V,(0l  ^ 


(B.ll) 


d_T 

dx 


(B.16) 


Ihe  derivative  dT(x),'dx  can  now  be 
evaluated  quite  easily.  Expanding  dT/dx 
in  terms  of  dT/dOj)  and  dx/dCiQ,  and 
using  equations  (B.9)  and  (B.11)  yields 

,T  dT  d 
dx  dx/d 


sm- 

v'M 


(B.  12) 


Substituting  equation  (B.16)  into  (B.2b) 
yields 


(6.17) 


The  limit  as  x  0  in  equation  (B.17)  is 
identically  the  same  as  the  limit  as 
Oq  ■*  0.  Incorporating  this  identity  and 
the  moments  definition  (see  equation 
(B.14)|  results  in 


III.  Determining  Cq 

(oeffitient  Cq  is  simply  the  normal 
incident  reflection  path  (source  to 


M 

c  =  •  (B.18) 

V.  Determining 

From  equation  (B.2c),  coefficient  C4  is 
given  by 

^  dx'  \dx' / 


To  evaluate  d ( dT^/ dx^ ) / dx^ ,  let  R2 
be  defined  as  dT^/dx^.  This  allows 
d ( dT^/dx^ ) / dx^  to  be  expressed  by 


(B.20a) 


d  /dT'\ 


dx‘  \dx 


2V2(0) 


■1.1 


X  dx/d'j. 


Substituting  for  x  and  dx/d6o  from 
equations  (B.IO)  and  (B.11)  and  re¬ 
arranging  results  in  the  following 
expression  for  d(dT2/dx2) /dx^ 


Taking  the  limit  as  x  goes  to  zero  and 
substituting  into  equation  (B.2c) 
results  in  the  following  expression  for 
C4  in  terms  of  moments. 


Expanding  dR2/dx  in  terms  of  dR2/d0o 
and  dx/d0o  al I ows 


j  dR,/dOp 
2  X  dx/dujj 


(B.20b) 


The  variable  R2  is  defined  as  dT^/dx^ 
which  is  related  to  "Ej  notation  in 
eouation  (B.16).  Incorporating  the 
"Ejj"  notation  all ows 


{B.20c) 


Expand i ng 
f  0  I  lowing 


equation 


results 


in  the 


jU _ 

x  dx/d'ip 


4M 


1  L 


1  - 


(B.22) 


VI.  Determining 

From  equation  (B.2d),  coefficient  C5 
is  given  by 


Cf. 


Lim  1  d  /  d  dT'  \ 
^  dx‘  ydx'  dx"  j 


{B.23) 


Defining  R4  as  d(dT2/dx2)/dx2, 
which  is  given  by  equation  (R.21  )  , 
allows  the  derivatives  of  equation 
(B.23)  to  be  expressed  as 


Following  the  s ame  steps  as  in  the 
previous  section  from  equations  (B.20a) 
to  (R.20c),  results  in  the  follov/ing 


d 

where  superscript  '  denotes  d6o  . 
Substituting  for  Eli  1  1 

(see  equation  (B.4))  results  in 
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Rearranging  anr)  substituting  for  x  and 
dx/dOg  from  equations  (B.IO)  and  (B.ll) 
results  in 


dx^’  \,dx‘  dxV  ^  E,,,  E,^3 

•  [dt(E;;.  -  E-,.,  )](»■«■=) 


Performing  the  indicated  differentia¬ 
tion  with  respect  to  Sg  yields  the 
following  equation  for  d [ d ( dT^ /dx^) / 
dx^ ] / dx^ 


‘'2  "  M, 


,  M  M 

C  =  1  -  _ i 

^4  2  ^  ? 

4M,  Mj 


_L_ 

SM,*^  L  ^ 


-  Mj  M.^  -  M_ 


d  /  d  dT  _ 3^_ 

dx'  ydx’  dx'  /  4  E  j  I  E'^ 


E=  -  E  , 


-El  ,  E, 

*  »  •  -  * 


E,  .  E  ,  ,  E, 


l.i  +  ^5  -1,1  T',,, 


(B.25) 

Taking  the  limit  as  x  anproachs  zero, 
and  substituting  into  equation  (R.23) 
yields  the  following  equation  for  C5 


Mj  +  24M3 


24M3’  M_j 


-  4M^  Mj’  -  5M^ 


M.i  M, 


where  the  moments,  M: ,  are  defined  by 


^  -I  . 


u  3^ 


1  [EM  M  , 

- - -  -  M.  M. 

sm/^L  ‘  ’ 


(B.26) 


VII.  Determining  C0 


C; efficient  Cg  can  be  derived  using  a 
procedure  identical  to  that  for  lower 
order  coefficients.  For  brevity,  only 
the  result  is  shown  here 


----  9M,’m 

5M,  L  ‘ 


,  +  24M,  M_| 


24M,  M 


4M,  M  '  -  5M 


(B.27) 


VIII.  Summary 

The  polynomial  moveout  coefficients  are 
related  to  moments  by  the  follov/ing 
express  ions 
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